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The one- dimensional electron gas exhibits spin-charge separation and power-law spectral responses 
to many experimentally relevant probes. Ordering in a quasi one-dimensional system is necessarily 
associated with a dimensional crossover, at which sharp quasiparticle peaks, with small spectral 
weight, emerge from the incoherent background. Using methods of Abelian bosonization, we derive 
asymptotically correct expressions for the spectral changes induced by this crossover. Comparison 
is made with experiments on the high temperature superconductors, which are electronically quasi 
one-dimensional on a local scale. 



In this paper, we consider the spectral signatures of 
dimensional crossover in the continuum theory of a quasi 
one-dimensional superconductor. This problem is of in- 
terest in its own right and for application to materials 
which are structurally quasi one-dimensional, such as the 
Bechgaard salts (organic superconductors). We believe 
that it is also interesting as a contribution to the the- 
ory of the high temperature superconductors. Although 
structurally these materials are quasi two-dimensk>nal, 
there is both theoretical and experimental evidenceEJ of a 
substantial range of temperatures in which "stripe" cor- 
relations make the electronic structure locally quasi one- 
dimensional, a phenomenon we have labeled "dynamical 
dimension reduction." Similarly, the ET versions of or- 
ganic superconductors are two-dimensional doped anti- 
ferromagnets, which we expect to show similar behav- 
ior. More generally, the high temperature superconduct- 
ing state emerges from a non-Fermi liquid normal state, 
often with a normal-state pseudogap. The quasi one- 
dimensional superconductor is the only solvable case in 
which such an evolution can be traced, theoretically. 

A quasi one-dimensional system can be thought of as 
an array of "chains" , in which the electron dynamics 
within a chain is characterized by energy scales large 
compared to the electronic couplings between chains. 
Since a one-dimensional system cannot undergo a finite 
temperature phase transition, any ordering transition 
with a finite critical temperature T c is necessarily associ- 
ated with a dimensional crossover. The electronic prop- 
erties at temperatures (or energies) large compared to T c 
can be understood by ignoring the interchain coupling, 
while at lower temperatures or energies, the behavior is 
that of a three dimensional system. 

In the one dimensional electron gasi (1DEG), as a 
consequence of spin-charge separation, the elementary 
excitations are collective modes with unusual quantum 



numbers and topological properties: The charge excita- 
tions are best understood as sound-like density-wave pha- 
sons (or, in dual representation, superconducting quasi- 
Goldstone modes) when the system is gapless, and as 
charge solitons with charge ±1 and spin when a charge 
gap is induced. (The precise meaning of the soliton 
"charge" is a quantized unit of chirality; See Eq. ( |33| ) 
and the subsequent discussion.) Similarly, the spin ex- 
citations of a spin-gapped system are spin solitons with 
charge and spin 1/2. When the elementary excitations 
do not have the quantum numbers of the experimentally 
accessible excited states, spectral functions do not exhibit 
sharp peaks corresponding to a well defined mode with 
a definite dispersion relation, lu = e(k). The single hole 
spectral function, G < (k,uj), which is measured in angle 
resolved photoemission spectroscopy (ARPES), involves 
excited states with charge e and spin 1/2, which thus 
consist of at least one charge soliton and one spin soliton. 
The dynamic spin structure factor S(k,u>), measured by 
neutron scattering, involves excited states with spin 1, 
which thus consist of two spin solitons. (We will see 
in Section [VC that, in fact, the relevant excited states 



contain two spin solitons and at least two charge antisoli- 
tons.) 

Below T c , where the system is three dimensional, we 
will show that the solitonic excitations of the 1DEG are 
confined in multiplets with quantum numbers that are 
simply related to those of the electron. For the case of 
three dimensional charge-density wave ordering, this has 
been known for some time. For the case of the supercon- 
ductor, it is related to the fact, noted recently by Salkola 
and Schriefferjj that either a spin soliton or a charge 
soliton induces a 7r-kink in the superconducting correla- 
tions. As a consequence of confinement, there is a finite 
probability of creating a final state consisting of a single 
bound spin and charge soliton pair in an ARPES exper- 
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iment. This will show up as a coherent (delta- function) 
piece in the zero temperature G < (fc,w). 

In this paper, we show that the coherent piece of the 
single particle spectral functions has a weight which van- 
ishes in the neighborhood of T c in proportion to a positive 
power of the interchain Josephson energy. It is this fact, 
that the spectral weight of the coherent piece is strongly 
temperature dependent below T c , rather than either the 
energy or the lifetime of the normal mode, which is a 
new feature that emerges from our analysis. Lh-js highly 
reminiscent of behavior observed in ARPESuu and in- 
clastic neutron scatteringQ measurements on the high 
temperature superconductors. We have also identified 
a new resonant feature in the spin spectrum of a quasi 
one-dimensional superconductor that emerges at temper- 
atures well below T c . 

If the 1DEG remains gapless down to T c , the super- 
conducting transition is BCS-like, in the sense that both 
pairing and phase coherence occur at the same time. In 
this case both are induced by the interchain Josephson 
tunnelling. We will mainly be concerned with the case 
in which a sort of "pairing" , i. e. the opening of a spin 
gap A s > 0, occurs in the Id regime well above T c . In 
this case T c is primarily associated with phase-jsrdering, 
and its scale is set by the superfluid density,0i3 rather 
than by the zero temperature single particle gap scale, 
Ao/2. In such circumstances the superconducting state, 
even at very low temperatures, maintains a memory of 
the separation of charge and spin which is a feature of 
the (Id) normal state. The unique "coherence length" of 
a BCS superconductor is replaced by two distinct corre- 
lation lengthsa: a spin length, £ s = v s /A s , where v s is 
the spin velocity, and a charge length, £ c — v c /A c , where 
Ac ~ 2T C . 

The remainder of the paper is divided into two self- 
contained parts; in Sections [j IV we derive asymptoti- 
cally exact results for the spectral properties of a quasi 
one-dimensional superconductor in the limit of weak in- 
terchain coupling. In Section [v| we summarize the prin- 
cipal results and discuss their application to experiment, 
especially in the high temperature superconductors. The 
reader who is interested only in results, not their deriva- 
tion, can skip the intervening sections. 

The model we study is defined by the Hamiltonian 



H 



E 
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dxHj + Hj 



(1) 



Hj is the Hamiltonian of 



where the sum runs over chains 
the 1DEG on chain j, and Hj is the Josephson coupling 
between chains. In Sees. [j III, we consider the single 
chain problem (Hj — 0). The problem is formulated 
using Abelian bosonization in Sec. Q. Next, we discuss 
the spectral functions for the 1DEG without (Section ||) 
and with (Section III) a spin gap - explicit expressions 
for various quantities in the presence of a spin gap are 
reported here for the first time of which we are aware. In 
Section IV, we extend these results to the case in which 



the most relevant interchain coupling is the Josephson 
tunnelling. An adiabatic approximation, which is exact 
in the limit where A s ^> A c , replaces the spin-charge sep- 
aration of the purely Id problem as the central feature 
of the spectrum - this section contains our principal new 
results. Applications to high temperature superconduc- 
tors are described in Sec. M. Various appendices expand 
upon the derivations in Section [V. 



I. ABELIAN BOSONIZATION AND THE 
SPECTRAL FUNCTIONS 

We begin by considering the properties of a single chain 
in the absence of any interchain coupling; we treat this 
problem using Abelian bosonization, which is based on 
the fact that the properties of an interacting 1DEG at low 
energies and long wavelength are asymptotically equal to 
those of a set of two independent bosonic fields, one rep- 
resenting the charge and the other the spin degrees of 
freedom in the system-. The widely discussed separation 
of charge and spirJl3'L-l in this problem is formally the 
statement that the Hamiltonian density Hj can be ex- 
pressed as 



H = H c + H S 



(2) 



where the chain index is implicit, and the charge and spin 
pieces of the Hamiltonian are each of the sine Gordon va- 
riety, 



H a — 



K a {d x 6 a f 



(d x (j) a ) 



- 1 



K a 



V a cos(v87n/> Q ) 



(3) 



where a = c, s for the charge and spin fields, respec- 
tively, 9 a is the dual field to <f> a , or equivalently d x Q a 
is the momentum conjugate to <p a - We consider a suffi- 
ciently incommensurate 1DEG and therefor set V c = 
since it arises from Umklapp scattering. Of course, if the 
Umklapp scattering is crucial to explain doped insulator 
behavior, its role cannot be neglected. Where there is no 
spin gap, or at temperatures large compared to A s , we 
can likewise set V s = 0. 

When V s is relevant, (perturbatively, this means K s < 
1) the spin gap is dynamically generated, i.e. it de- 
pends both on V s and the ultraviolet cutoff in the 
problem, A, according to the scaling relation A s ~ 
v s A[V s /v s A 2 ] 1 ^ 2 ^ 2Ks \ At the gapless fixed point, spin- 
rotational invariancc requires K s = 1, at which point V s 
is perturbatively marginal. It is marginally irrelevant for 
repulsive interactions (K s > 1) and marginally relevant 
for attractive interactions (K s < 1). Thus, the long dis- 
tance spin physics is described by H s with V s = and 
K s = 1 for a gapless spin-rotationally invariant phase. 
Where there is a spin gap in a spin rotationally invariant 
system, it is exponentially small for weak interactions, 
A s ~ iJV s v s exp[-v s A 2 /2ttV s }. 
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In order to compute correlation functions, we use the 
Mandelstam representationlij of the fermion field opera- 
tors 

ip\,o-(x) — No exp [iXkp x - i$ A)(r (x)] , (4) 

where M a contains both a normalization factor (which 
depends on the ultraviolet cutoff) and a "Klein" factor 
(which can be implemented in many ways) so that Af a 
anticommutes with N a > for a ^ a 1 and commutes with it 
for a = a' . In addition 

*A, ff = [(Pe - He) + v{e s - \<t>s)\ , (5) 

where A = — 1 for left moving electrons, A = +1 for 
right moving electrons, and a = ±1 refers to spin po- 
larization. From Eq. (||), it is a straightforward (and 
standards) exercise to obtain the boson representations 
of all interesting electron bilinear and quartic operators. 
Physically, <f> c and <f> s are, respectively, the phases of the 
2kp CDW and SDW fluctuations, and 9 C is the super- 
conducting phase. The long-wavelength components of 
the charge (p) and spin (S z ) densities are given by 

^— ' ' n \ it 

A, (7 

A,er 7F 

When analyzing results for this model, it is always im- 
portant to remember that the parameters which enter the 
field theory are renormalizcd, and are related to the mi- 
croscopic interactions in a very complicated manner. For 
instance, although for a single-component 1DEG with 
repulsive interactions V s is always irrelevant, for mul- 
ticomponent IDEG's, and for the "1DEG in an active 
environment" , it is common to find a dynamically gen- 
erated spin gap, even wiiea-the microscopic interactions 
are uniformly repulsive .E3~li3 

The bosonized expressions for all electron operators 
are readily extended to an array of chains by adding a 
chain index to the Bose fields and to the Klein factors; 
the Klein factors on different chains must now anticom- 
mute with each other. Where single particle interchain 
hopping is relevant, the Klein factors appear explicitly in 
the bosonized Hamiltonian. Where only pair hopping 
and collective interactions between neighboring chains 
need be included in the low energy physics, the Klein 
factors cancel in H. 

While it is generally simpler to derive results concern- 
ing the spectrum, it is important for comparison with ex- 
periment to compute actual correlation functions. Specif- 
ically, we will consider the transverse spin dynamic struc- 
ture factor 

S(x,t;T) = (S x 2kp Hx,t)S x 2kF (0,0)) 

+ <Sf fcF t (x,i)S| fcF (0,0)), (6) 



where 

(7, a' 

and the r are Pauli matrices. We will also consider the 
one-hole Green function, 

G<(x,t) = (^_^(x,t)^ 1A (0,0)) , (8) 
the singlet-pair correlator, 

x(x,t) ee (4 t (m)VA u (M)V--u(o,o)^t(o,o)) , (9) 

and the various spectral functions, S, G < , and Xj ob- 
tained by Fourier transforming these correlators. As a 
consequence of separation of charge and spin, S, G < , 
and x are expressible as a product of spin and charge 
contributions, and consequently 5, G < , and x are con- 
volutions. For instance, 

G<(k,co) = J ^^- Gs {k-q,u:-v)G c {q,v). (10) 



II. HIGH TEMPERATURE: LUTTINGER 
LIQUID BEHAVIOR 

At temperatures large compared to T c and the spin 
gap, A s , (or at all temperatures in systems in which 
T c = A s = 0), the 1DEG exhibits "Luttinger liquid" 
behavior. Because the Luttinger liquid is a quantum crit- 
ical system, the response functions have a scaling form. 
Specifically, this impliescl that 

G<(k,uj;T) = T 2 ^ +2 ~<°- 1 G < (k/T,uj/T:l) , (11) 

where we define for a = c or s 

l a = \{K a + K- x -2), (12) 
and that so long as K s > 1, 

S(k, lo; T) = T^ 1+K ^S{k/T, u/T; 1) . (13) 

Note that here, and henceforth, we will measure k relative 
to kp and 2kp, respectively, when computing the scaling 
functions G < and S. If the system is spin-rotationally 
invariant, K s = 1 in the above expressions. 

The form of these scaling functions can be computed 
analytically in many cases; this has recently been accom- 
plished in Ref. ^0|. They may or may not have a peak 
at energies small compared to the bandwidth, depending 
on certain exponent inequalities. Where there is a peak, 
it occurs at positive energies u> = ±v a k + (const. )T, but 
the peak width, however defined, does not narrow in pro- 
portion to T at low temperatures; such a peak does not 
correspond to a quasiparticle. 
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III. INTERMEDIATE TEMPERATURE: THE 
LUTHER-EMERY LIQUID 

When V s is relevant, the spin sine Gordon theory scales 
to a strong-coupling fixed point, and the excitations are 
massive solitons, in which <p s changes by ±-y/7r/2 (i.e. 
S z = ±1/2). This problem is most simply treated in 
terms of spin fermion fields, 



(14) 



The refermionized form of the Hamiltonian is then 

Us = «;.[¥+ , - vl A d x v Sil } 

+A,[*J +H.C.] 



(15) 



where 



Vs = 

A s = 
9s = 



4K S 



K, 



A 

27TU, 



4K X 



K, 



(16) 



Forp£C s = 1/2, which is known as the Luther-Emery 
point J13 the refermionized model is non-interacting and 
massive, with a gap A s = A s . Assuming there is 
a single massive phase of the sine Gordon theory, the 
Luther-Emery mode£3 will exhibit the same asymptotic 
behavior as any other model in this phase. Formally, 
the Luther-Emery point can be thought of as a strong- 
coupling fixed point Hamiltonian, and g a , which vanishes 
at the fixed point, is the amplitude of a leading irrelevant 
operator .Ell We will henceforth compute correlation func- 
tions at the Luther-Emery point, and then comment on 
the effects of deviations from this point. 

Now, in computing the various spectral properties of 
the system, we can distinguish two regimes of temper- 
ature: at temperatures large compared to A s , the spin 
gap is negligible, and the results for the Luttinger liq- 
uid apply. If the temperature is small compared to the 
spin gap, then we can compute the spin contributions to 
the various correlation functions in the zero-temperature 
limit, and only make exponentially small errors of order 
exp(— A 8 /T). The spin piece of the transverse spin re- 
sponse function can be expressed in terms of the spin 
fermion fields 



(17) 

Since the theory reduces, at the Luther-Emery point, to a 
theory of free massive fermions, the corresponding spec- 
tral function can be readily computed with the result, for 
T = 0, 



Avj\q 1 E s (q 2 ) - q 2 E s (q 1 )\ 
where the spin soliton spectrum is 



(18) 



E s (k) = ^vW+M 



(19) 



and q\ t 2 are the solutions of the quadratic equation 
uj + E s (q) + E s (k + q) = 0. Explicitly 



4A? 



v 2 k 2 — uj 2 



(20) 



The spin piece of the one hole Green function is more 
complicated, since it involves nonlocal operators in the 
refermionized form: 

G s (x,t) = (itffoijtftfot)*.,-! (0,0)^.(0,0)) , (21) 



where the vertex operator U s (x) 



W2<M*) with 



<f) s (x) 



/ 7r/2 Ylx J X dy^\ \^s,\- From kinematics, it 



follows that this Green function consists of a coherent one 
spin soliton piece and an incoherent multisoliton piece: 

G.(k,u) = Z s {k)5[u - E s (k)} + G s mum \k,u) , (22) 

where the multisoliton piece is proportional to Q[lu — 
3i? s (fc/3)]. (Deviations from the Luther-Emery point in 
the case g s > will result in the formation of a spin 
soliton-antisoliton bound state, a "breather" , which can 
shift the threshold energy for multisoliton excitations 
somewhat.) 

At the Luther-Emer*-, point it is possible to obtain 
closed form expressionsE3 for the matrix elements of the 
vertex operator between the vacuum and various multi- 
soliton states and from that to compute Z s explicitly. We 
will report this calculation in a forthcoming paper, Ref. 
pfj| . Here, we use a simple scaling argument, which can 
be generalized to the case of nonzero interchain coupling, 
to derive the principal features of this result, especially 
the dependence of Z s on A s . In the absence of a spin 
gap, and at T = 0, G s can be readily evaluated to give 
the scaling form 

7T (V S A)5- 2 T» ,vy,-l/ , w.-i 



r(7-)r(7. + h) 



x0(w- v s \k\) . (23) 



Because the sine-Gordon field theory is asymptotically 
free, the high energy spectrum, and hence the depen- 
dence of G s on A, is unaffected by the opening of a spin 
gap. With this observation, it is simply a matter of di- 
mensional analysis to see that 



z 8 (k) = (A£ s )*- 2 ^/ s (fc£ s ) 



(24) 
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where £ s = v s /A s is the spin correlation length. f s is a 
scaling function which is independent of K s . It can be 
calculatedc2l using the exact matrix elements available for 
K s = 1/2, with the result 



(25) 



f.(x)=c[l- 



where c is a numerical constant. r-. 

The abova extends the earlier results of Voitu and 
Wiegmann.Ej In particular, the analytic structure (as a 
function of k and tu) of the one soliton contribution to 
Eq. p3) reproduces that found in earlier work. The ex- 
plicit discussion of the nonanalyticities due at the three 
soliton threshold is new, although fairly obvious; more 
muted singularities occur at the five and higher multi- 
soliton thresholds, which we will not discuss explicitly. 
The specific expression in Eq. ( |24| ) is, to the best of our 
knowledge, new, and is the most important feature of 
this result for the purposes of the present paper. 

The charge pieces of both response functions are unaf- 
fected by the opening of the spin gap. Consequently, S 
and G < have power law features (which can be a peak or 
a shoulder depending on K c ) at u> = 2E s (k/2)+0(T) and 
to = E s (k)+0(T), respectively, with a shape and temper- 
ature dependence, both readily computed, determined by 
the still gapless charge density fluctuations. For exam-, 
pie, we can evaluate the spectral function explicitlyc2l 
at T = in the limit V s /v c — > (and for arbitrary 
uj < 3E s (k/3)), or when \u — A s \ <C A s (for arbitrary 
v s /v c ): 



Av, 

4r( 7c )r( 7c + i) [~2 

x Z s {k)[u-E s {k)] 2 ^ 



-3-27c 



*9[u) - E a {k)] . (26) 



Here B(x,y) is the beta function. Again, the fact that 
these excitations are not quasiparticles is reflected in the 
fact that, even where peaks in the spectral function oc- 
cur, they do not narrow indefinitely as T — > 0. 

In the presence of a spin gap, the spin contribution to 
the long distance behavior of the superconducting sus- 
ceptibility is a constant; 

X(x,t) ~ \(u^\ 2 (e iV ^ 6 ^e- iV ^ 6 ^) 

~ (A£ s )-Ks( e i^e c (x,t) e -iV2^e c (o,o)} . (27) 

From this, one sees that, within a chain, one can identify 
(A£ s y Ks/2 as the "amplitude" and \/2tt9 c as the "phase" 
of the order parameter .113 



IV. LOW TEMPERATURE: THE 3D 
SUPERCONDUCTING STATE 

For temperatures of order T c and below, interchain 
couplings cannot be ignored. Single particle hopping and 



all magnetic couplings are irrelevant by virtue of the pre- 
existing spin gap. For K c > 1/2, the Josephson coupling 
is perturbatively relevant, but for K c < 1, the 2kp CDW 
coupling is more relevant. For the simplest realizations of 
the 1DEG, K c < 1 corresponds to repulsive interactions 
between charges. However, we have showrc£lli3 that for 
fluctuating or meandering stripes, such as occur in the 
high temperature superconductors, the CDW coupling 
gets dephased, so that the Josephson coupling is the most 
relevant, even when 1/2 < K c < 1. 

Since we are interested in the onset of superconductiv- 
ity, we consider the case in which the Josephson coupling 
between chains is more relevant. The pair tunnelling in- 
teraction between chains, which appeared in Eq. (0), can 
be simply bosonized: 



Hj = -J sc J2 J ^[AtAj + H.G 



(28) 



<i,3> 



where the pair creation operator on chain number j is 

A*(x,t) = Vi, T ^I u +V>li, T ^u 

oc cos(-\/27r(^ s ) exp(iv / 27T0 c ) , (29) 

and we have left the chain index implicit. 

Since the state below T c has long range order, and since 
we assume that the coupling between chains is weak, it 
is reasonable to treat it in mean field approximation^ 
although we continue to treat the one-dimensional fluc- 
tuations exactly. Thus, rather than considering a full 
three-dimensional problem, we consider the effective sin- 
gle chain problem defined by the Hamiltonian 



H = H S +Hc-J cos(V27r<j) s ) cos(V2tt6> c ) , (30) 

where J is related to the pair tunnelling amplitude by 
the mean field relation 



(31) 



J = zJ sc (A/7r) 2 (cos(V27T0 s )cos(V27r<9 c )) 



where z is the number of nearest neighbor chains. [Since 
the average of cos(v / 2~7r0 s ) sin(v / 27T0 c ) vanishes, no sine 
term appears in the effective Hamiltonian (|30|).] Note 
that the pair hopping term in Eq. (|3^) couples charge 
and spin, as is characteristic of higher dimensional cou- 
plings. 

The mean field approximation is exact in the limit of 
large z and small zJsc- In three dimensions, this mean 
field approximation will produce some errors in the crit- 
ical regime in the vicinity of T c , but because of the long 
correlation length along the chain just above T c , the crit- 
ical region is always small for small Jsc: and well below 
T c , this approximation is safc£3 

Because of the presence of relevant cosine terms, there 
are superselection rules which divide Hilbert space into 
various soliton sectors. The soliton sectors are specified 
by two integrals: 
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N s = v/2/tt / dxd^s (32) 

J — oo 

= [<f> 3 (oo) - ^ (-00)] = 2 / efeS^ , 



and 



A c = v/27^ / f£c<9 x c = ^2/^ [6> c (oo) - c (-oo) 



(33) 

N s is simply the number of spin solitons minus the 
number of antisolitons or the total value of S z in units of 
h/2. The interpretation of N c is a bit more subtle. Since 
we are looking at a superconducting state, the r&lectro- 
static charge of a quasiparticle is not defined.EStiij How- 
ever, N c is a conserved "chirality" equal to the number of 
right moving minus the number of left moving electrons, 
so that we can still interpret eN c as a sort of quasiparticle 
"charge" ; it represents .the coupling of the quasiparticles 
to a magnetic fluxOEao 

The presence of the cos(v / 87n/>s) term in the single 
chain Hamiltonian results in the quantization of N a in in- 
teger units. The presence of the cos(^/2tt(/) s ) cos(v / 2~7r0 c ) 
term in Ti results in the quantization condition that 
N s + N c be an even integer! Physically, this means that 
excitations can have spin h and charge (N s = 2 and 
N c = 0), spin and charge 2 (N s = and N c = 2), 
spin h/2 and charge 1 (N s — 1 and N c = 1), etc., but 
that all the exotic quantum numbers of the soliton ex- 
citations of the isolated 1DEG are killed. Formally, the 
addition of the pair hopping term to the Hamiltonian of 
the 1DEG leads to a confinement phenomenon. Along 
the entire segment of chain between two spatially sepa- 
rated ±y/n/2 solitons, there is a change in sign of the 
pair hopping term. (See Eq. (|29]).) This leads to an en- 
ergy which grows linearly with the separation x between 
solitons, ~ J\x\, regardless of whether they are charge 
or spin solitons or antisolitons. 

The importance of this observation becomes clear when 
we study the operators in whose correlation functions we 
are interested. Since 



My) - 



V2 9(y-x) , 
(34) 



and 



--e c {y) + ^/2Q{x-y) , 
(35) 

it is clear that the fermion annihilation operator ^— i,f 
creates a spin antisoliton and a charge antisoliton, while 
the 2kp piece of the spin-raising operator, 5 , ^ Cf , creates 
a pair of spin solitons and a pair of charge antisolitons. 
Both these combinations decay into a set of free solitons 
in the absence of the interchain coupling, but in its pres- 
ence, the former becomes a bound state, and the latter 



a resonant state. Thus, G < develops a coherent piece 
with a well defined dispersion relation as superconduct- 
ing phase coherence between chains occurs. S develops 
a resonant peak at a temperature well below T c . 



A. Zero spin soliton sector 

For the case in which the spin gap A s of the isolated 
chain is large compared to the interchain coupling, the 
fluctuations of the spin field are high energy (fast) com- 
pared to any charge fluctuations, and indeed only slightly 
affected by the onset of superconducting order. In this 
limit, the eigenstates can be treated in the adiabatic ap- 
proximation. 

In the ground state (N s = 0) sector, the spin field fluc- 
tuations are little affected by Hj; all spin correlations 
can thus be computed as in the previous section. More- 
over, because of the spin gap, so long as T -C A s , the 
spin fields can be approximated by their ground state. 
For computing the charge part of the wave function, we 
can replace the operator cos(v / 2~7r</ ) s) in Hj by its expec- 
tation value at zero temperature in the decoupled ground 
state, 



cos(V27r<^) -» (cos(V27r<£ s )) = C s - (A£ s )" 



(36) 



where the subscript "o" refers to the expectation value 
in the ensemble with J$c set equal to zero, (see also Ref. 
|32] ). This leaves us with a sine Gordon equation for the 
charge degrees of freedom, with potential 



(37) 
(38) 



JC S cos(V2tt0 c ) . 
Again, we solve this problem by refermionizing 



:A/" c exp[zv/V2(6» c -2A(/> c )] 



The refermionized form of the Hamiltonian is 

n c = «c[*t-A*=.-i ~ *lA*c,i] 

-A c [*t *t +h.c] 



+flc*c,l*l ,_i*c,-l*c,l 



(39) 



where 



1 



K, 



ttJCs 
A 



g c = 2irv c [ — - K L 



(40) 



Since N s = 0, the superselection rule implies N c = 2m, 
which upon refermionization is simply the condition: 



]TA f dx[^\^ c , x }=N c /2 



(41) 



G 



It is also interesting to note that the superconducting 
pair creation operator can be expressed in an intuitively 
appealing form in terms of charge soliton creation oper- 
ators 



(42) 



Recall that here the charge solitons are spinless fermions. 
This expression emphasizescj the fact that spin gap for- 
mation, which is associated with the quenching of the 
fluctuations of the spin density phase, (f> s , can also be 
identified with the growth of the amplitude of the super- 
conducting order parameter. While the charge solitons 
clearly also make a contribution to the amplitude of the 
order parameter, the phase of the order parameter comes 
entirely from the charge. 

For K c — 1/2, just as for the Luther-Emery point 
for the spin fields, the refermionized Hamiltonian for 
the charged excitations is noninteracting and massive 
(gapped), and A c = A c . In computing the asymptotic 
form of correlations we will set K c = 1/2. We can now 
readily compute the expectation value of the pair hop- 
ping term so as to relate two physically important quan- 
tities: the excitation energy scale, A c , and the interchain 
portion of the internal energy 



A r 



\ C,l C, 



H.c.)= l 7(cos(v / 2~7r0 s ) cos( 
= (A c /7r£ c )u (A c ,T) 



2tt0 c )) 
(43) 



wherefc = v c /A c is the charge correlation length. Equa- 
tion (|4^) has the form of a BCS gap equation with 



pv c A 

u (A c , T) — I dx 
Jo 



1 



tanhl —y/ x 2 + A 2 c , 

(44) 



where the mean field relation for A C (T) is 

TTV r 



u (A c ,T) 



zJscC 2 s ' 

Consequently we find the familiar BCS relations 
T c = 0.57A C (0) , 

A c (0) = 2v c A exp[-nv c /zJ SC C 2 s } , 



(45) 



(46) 



(47) 



As is shown in Appendix B, for a system with equal 
areas of domains in which the stripes run along the 
x and y directions, the macroscopic phase stiffness is 
equal to the geometric mean of the superfluid density in 
the directions parallel and perpendicular to the chains, 
k = y^Kjj k±. Since the phase stiffness along the chains 
is simply K|| = v c K c , it follows that R(T = 0) is (up to 
logarithmic corrections coming from uq) simply propor- 
tional to T c . This is a microscopic realization of a more 
general phenomenon which occurs in systems with low 
superfluid density p it is phase ordering, as opposed _to 
pairing, which determines T c . In a future publications, 
we will study the effects of quantum and thermal phase 
fluctuations on the evolution of the superfluid density of 
a quasi one-dimensional superconductor. 

With little additional effort, we can study the pair field 
susceptibility \ a t energies small compared to 2 A,,. In 
this low energy limit, as in Eq. (p7|), we can replace 
the spin operators in \ by their ground state expectation 
values. 

The charge part of x can be expressed in terms of the 
charge fermion fields 



(50) 



At the free charge fermion point (K c — 1/2) the corre- 
sponding spectral function is readily evaluated, for T = 
and oj <C 2A S , with the result 



+ 



C s u 



4£ c 2 (fc/2) 



2A? 



Av 2 \ qi E c (q2) - q 2 E c (q 1 )\ 



Q[uj -2E c (k/2)] 



(51) 



where E c (k) and g 12 are the analogs of Eqs. (19 2^) with 
A c substituted for A s and v c for v s . 

Away from the Luther-Emery point, if g c > (K c < 
1/2), the two solitons repel, and hence the effect of g c 
can be ignored, but for g c < (K c > 1/2), there is an at- 
tractive interaction between the two solitons and hence, 
this being after all a one dimensional problem, they form 
a bound state. This will slightly modify the expression 
for x- 



A C (T) « 1.74A c (0)-\/l -T/T c for T « T c . (48) 

In general, the actual form of A c (0) in terms of Jsc an d 
C s is modified according to the microscopic value of K c . 

The transverse superconducting phase stiffness n± 
(proportional to the superfluid density) is 

kl = 2na(Hj) , (49) 

where d is the spacing between chains and (Hj) is given 
in Eq. Thus at zero temperature k± ~ T 2 /v c . 



B. The one hole sector 

In the one soliton sector of the spin Hamiltonian, the 
adiabatic approximation requires reexamination. While 
for the most part, the spin modes are fast compared to 
the charge modes, the Goldstone mode (translation mode 
of the spin soliton) is slow compared to all other modes, 
and so must be treated in the inverse adiabatic approxi- 
mation. Thus, we consider the charge Hamiltonian with 



7 



a spin antisoliton at fixed position, R s . The pair tun- 
nelling term is then 



(52) 



JC s s\gT\{x — R a ) cos(v27r6', 



where we have used the fact that £ s = v s /A s (which char- 
acterizes the width of the spin soliton) is small compared 
to the charge correlation length, £ c = v c / A c , to approx- 
imate the profile of the spin soliton by a step function. 
Upon refermionization, the charge Hamiltonian is still of 
the same form as Eq. (|9|) with the term proportional to 
A c replaced by 

A c ^ -A c sign(x-i? s ) . (53) 

For K c = 1/2, upon the canonical transformation, 

V>c,-i = , ^c,i = *c,i , (54) 



the charge soliton Hamiltonian is of the same fornO 
as the fermionic Hamiltonian of a commensurability 2 
Peierls insulator, "polyacetylene" , in-fhe presence of a 
topological soliton. As is well known,Ej there is an index 
theorem that implies the existence of a zero energy bound 
state associated with the soliton, the famous "midgap 
state" or "zero mode" . All other fermionic states have en- 
ergies greater than or equal to A c . Importantly, since in 
this sector N s = —1, the superselection rule N c — 2m+l, 
requires that the fcrmion number is half integer! 

- I dx : [^ a Vo,a] : = N c /2 = m + 1/2 . (55) 

This is essential, since_with the midgap state occupied 
the fermion number isE-3 +1/2, while with it empty the 
fermion number is —1/2. The midgap state is associated 
with the bound state of the spin and charge antisolitons. 

To compute the charge contribution to the soliton cre- 
ation energy we need to evaluate the difference between 
the ground state energies of the charge Hamiltonian in 
the presence and absence of a kink. We have done this 
by taking the limit of vanishing soliton width of_a gen- 
eral expression of Takayama, Lin-Liu, and Maki£j (and 
dividing by 2 for the spinless case) . The resulting soliton 
creation energy is just A c /2 ; in other words, the rest en- 
ergy of the electron, i. e. the bound state of a spin soliton 
and a charge soliton, is 



A 



A c /2 m A 5 



(56) 



From this discussion, we can immediately conclude 
that for T <C T c <C A s , the one hole spectral function 
has a coherent piece and a multiparticle incoherent piece, 

G < {k,w) = Z(k)8[w - 8{k)\ + G {mum) (k,Lu) , (57) 

where 



This follows from the fact that the bound state of a spin 
soliton and a charge soliton has the same quantum num- 
bers as a hole. The multiparticle piececS has a threshold 
slightly above the single hole threshold at ui = £(fc)+2A c . 

The overlap factor, Z(k), contains factors from both 
the spin and the charge parts of the wavefunction; so 
long as k£ s « 1, Z(k) = Z c (k)Z s (Q) where ZJO) de- 
pends on the spin correlation length as in Eq. (|24j) , and 
Z c {k) contains all remaining contributions. We can ob- 
tain a scaling form for Z c using the same method of anal- 
ysis employed previously for Z s . Specifically, at T = 
in the absence of interchain coupling, and for lu -C 3A s 
and |fc£ s | -C 1, G < is given by the expression in Eq. (gg). 
Since the opening of a charge gap does not affect the high 
energy physics, the dependence of G < on A is unaffected 
by the interchain coupling. Indeed, so long as A c <C A s , 
the dependence of G < on A s is likewise unchanged. Thus, 
by dimensional analysis, it follows that 



Z(k) = Z s (0)(A£ c )-i- 2 ^A 7c />£ c ) 
where / is a scaling function and 



A, 



g( 7c , 7c + l/2) 

r( 7c )r( 7c + 1/2) 



(59) 



(60) 



£(k) = 



A§ 



(58) 



Unfortunately, we do not have exact results from which 
to compute f{x) explicitly, but there is no reason to ex- 
pect it to have any very interesting behavior for small 
x. 

At temperatures between T = and T — T c , the same 
arguments lead to a simple approximate expression for 
the spectral function. Specifically, the principal temper- 
ature dependence comes from A c which is a decreasing 
function of T. At mean field level, the temperature de- 
pendence of A c can be computed from Eq. (Q) . In par- 
ticular it vanishes at T c according to Eq. ( |48| ) . Since fluc- 
tuation effects produce superconducting correlations be- 
tween neighboring chains at temperatures above T c , this 
simple mean field behavior will be somewhat rounded, 
but the qualitative point that A c becomes small at tem- 
peratures above T c is quite robust. 

Consequently, the quasiparticle weight, Z, which is 

proportional to A c c 2 , is a strongly decreasing func- 
tion of T which vanishes in the neighborhood of T c . 
The quasiparticle gap, A , on the other hand, is only 
weakly temperature dependent, dropping from its max- 
imum value A = A s + |A C (0) at T = to A = A s 
in the neighborhood of T c . Scattering off thermal ex- 
citations will, of course, induce a finite lifetime for the 
quasiparticle at finite temperatures. 

Neither a charge soliton nor a spin soliton can hop from 
one chain to the next, but a hole can. The problem of the 
transverse dispersion of the coherent peak in the single 
hole spectral function is addressed in Appendix A. Not 
surprisingly, we find that the effective interchain hopping 
matrix element, t±, is replaced by an effective interchain 
hopping matrix clement, 
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tf s = Z(k)t± 



(61) 



V(R) = A c [l + v(R/t[ c )], (67) 



Thus, the dispersion of the coherent peak transverse to 
the chain direction is an independent measure of the de- 
gree of interchain coherence. 

C. The two spin soliton sector 

To compute S, we need to study states in the N s = 2 
sector. Interestingly, (in contrast to the case of an or- 
dered CDW) in a quasi one-dimensional superconduc- 
tor, the 2fcp spin density wave operator also creates two 
charge antisolitons: iV c = —2. Again, for the most part, 
the spin fluctuations are fast and high energy compared 
to the scale of the charge fluctuations, and can thus be 
treated in the adiabatic approximation - indeed, they are 
little affected by the presence of the interchain Josephson 
coupling. However, there are two low frequency modes 
associated with the soliton translational degrees of free- 
dom, which must be treated in the antiadiabatic approxi- 
mation. Consequently, we obtain an effective Schrodinger 
equation governing the center of mass motion of the two 
spin solitons: 

H^^2A s -^±^ + V( Xl -x 2 ), (62) 

3=1 ■? 

where Xj is the position of soliton j, 

AP = A s /v 2 s , (63) 

and V is the adiabatic spin soliton potential, obtained 
by integrating out the (relatively fast) fluctuations of the 
charge degrees of freedom. 

To compute V(R), we again rely on the analogy be- 
tween the refermionized version of the charee.nart of the 
Hamiltonian and solitons in polyacetylene.tTE3 In this 
case, V(R) is recognized as the difference in the ground 
state energy of a massive Dirac fermion in the presence 
and absence of a pair of zero width solitons separated by 
a distance R, i.e the Hamiltonian in Eq. ( |39| ) with 

A c -» A c sign(4:r 2 - R 2 ) . (64) 

Since 2N C = —2, this energy difference is to be computed 
in the fermion number —1 sector. 

From the results in the previous section, it follows that 

V(R) -> A c as R -*■ oo , (65) 

since in this limit, the two solitons are noninteracting, 
and reduce to the solution discussed in the previous sec- 
tion. Similarly, since for R — 0, the energy approaches 
that of the uniform system with fermion number — 1, 

V(R) -> A c as i? -> . (66) 

Moreover, from simple scaling, it is clear that 



where v(x) is independent of the magnitude of A c and 
v(x) — > for x — > and x — > oo. For intermediate R/S, c , 
we have been unable to obtain an analytic expression for 
v, although it is easily derived numerically, as described 
in Appendix C, with the result shown in Fig. 1. As can 
be seen, v(x) rises from to a gentle maximum at x w 0.3 
where v(0.3) « 0.2, and then drops exponentially back 
to zero at large separation. 

What this means is that there is no true bound state 
in the spin one excitation spectrum. The spin one 
excitations, even in the superconducting state, are al- 
ways unstable to decay into a pair of far separated spin 
1/2 quasiparticles. However, near the threshold energy, 
lu = 2A S + A c , there is a nearly bound (resonant) state 
with a lifetime which is exponentially long. Treating Eq. 
( |62| ) in the WKB approximation, we see that the decay 
rate of the resonant state is 



r ~ exp [ - B{v c /v s )^/A s /A c ] , (68) 
where B is a constant of order 1. 
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FIG. 1. The scaled version of the adiabatic potential, v(x), 
defined in Eq. (67), computed from the SSH model in a sys- 
tem of size 3000 sites with open boundary conditions. The 
different curves are for different magnitudes of dimerization 
corresponding to a coherence length of the indicated magni- 
tude in units of the lattice constant. 

Using the fact that S ~ A~~ in the absence of a 
charge gap and utilizing the same scaling arguments ap- 
plied previously to the coherent piece of G < ) it is easy 
to see that the weight associated with this resonant state 

is proportional to A c Kc . However, because the barrier 
height is small compared to A c , the thermal decay of the 
resonant bound state will become large, due to activation 
over the barrier, at a temperature well below T c . 



<) 



D. The "BCS-like case": no preexisting spin gap 

When there is no spin gap on the isolated chain and 
there are repulsive interactions in the charge sector, the 
interchain Josephson coupling is perturbatively irrele- 
vant. Thus, the usual case for a quasi one-dimensional 
superconductor is the already analyzed case with a preex- 
isting spin gap. However, it is worthwhile considering the 
case (with K c < 1) in which both the spin gap and the su- 
perconducting coherence are induced by a relevant inter- 
chain Josephson coupling. This case, even though quasi 
one-dimensional, is much more akin to the usual BCS 
limit, in that there is a single gap scale in the problem, 
and pairing (gap formation) and superconducting coher- 
ence occur at the same temperature and with roughly the 
same energy scale. This .case has been analyzed exten- 
sively in the literatureJlJ Ej It should be noted, however, 
that here, too, since the "normal" state is a non-Fermi 
liquid, the coherent piece of all spectral functions will be 
strongly temperature dependent below T c , and vanish in 
the neighborhood of T c . 



V. SUMMARY AND IMPLICATIONS FOR 
EXPERIMENT 

In this paper, we have obtained explicit and detailed 
results for the properties of the superconducting state 
of a quasi one-dimensional superconductor. We have 
studied this problem as a quantum critical phenomenon, 
in which the quantum critical point is reached in the 
Id limit of no interchain coupling, and hence we have 
treated the interchain Josephson coupling as a small pa- 
rameter. In particular, we expect (as discussed below) 
that the results are pertinent to underdoped and opti- 
mally doped high temperature superconductors, where 
self-organized stripe structures render the system locally 
quasi one-dimensional. 

It is often argued that, even in fairly exotic circum- 
stances, and even when the normal state is a non-Fermi 
liquid, the superconducting state itself is fairly conven- 
tional and BCS-like. We have shown that there are a 
number of ways in which this expectation is violated. In 
the first place, there are two "gap" scales, A s ^> A c , 
whereas in a BCS superconductor there is one, Ao, and 
correspondingly two correlation lengths, £ s and £ c , in 
place of the one, £o, of BCS theory. However, both 
gaps are, in a very real sense, superconducting gaps: A s 
is associated with spin pairing (i.e. a nonzero value of 
(iff g i^s^-i)) and the existence of a local amplitude of the 
order parameter. A c is a measure of interchain phase co- 
herence. In the case in which there is a preexisting spin 
gap on the isolated chain, £ s remains finite at the quan- 
tum critical point, whereas £ c diverges. The same holds 
true in the superconducting phase, a bit away from the 
quantum critical point, where £ c diverges at T c , while £ s 
remains finite£a 



It is perhaps worth noting that many of these novel 
aspects of the superconducting state are considerably 
more general than the-.particular model we have solved. 
Indeed, recently, Leec3 derived similar results from the 
gauge theory formulation of a flux phase to superconduc- 
tor transition. While this derivation presupposes rather 
different seeming microscopic physics, it does build in the 
doped insulator character of the superconducting state, 
which is the essential feature of the results. Likewise, 
many features we have discussed here bear a close resem- 
blance to the dimensional crossover from a conjectured 
2d non-Fermi liquid to a 3d superconductor envisaged 
in the context of the inter-layer tunnelling mechanism of 
high temperature superconductivityo. 



A. Summary of Results 

For the benefit of the reader who skipped the technical 
exposition, we begin by summarizing our most important 
results. We consider here the case in which there is a pre- 
existing spin gap, A s /2 3> T c , on an isolated chain, and 
we focus on the effects of the interchain Josephson cou- 
pling between stripes at lower energies. 



1. Thermodynamic Effects 

The effect of the interchain Josephson coupling is to 
produce an interchain coherence scale, A C (T). At mean 
field level, A C (T) vanishes for any T above T c , and 
while fluctuation effects will produce a small amount of 
rounding to this behavior, because of the large coherence 
lengths along the chain the degree of rounding will al- 
ways be small in the quasi- ID limit. It is the coherence 
scale that determines T c , in the sense that 

T c » A c (0)/2 « A s /2 . (69) 

(A c is expressed in terms of the strength of the Josephson 
tunnelling matrix elements in Eq. ([47]).) The superfluid 
densities in the directions transverse and parallel to the 
chain direction are, respectively, 

k_l = 2au Al/v c , k\\ = v c K c , (70) 

where uq is a constant (see Eq. (|44|)) which depends 
weakly on parameters, d is the spacing between chains, 
v c is the charge velocity, and K c is the charge Luttingcr 
parameter. In two dimensions, if, on average, there is a 
4-fold rotationally invariant mixture of domains in which 
the chains run along the x and y directions, respectively, 
the macroscopic superfluid density is isotropic and given 
by 

R(T) = A C (T) . (71) 
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2. Single Hole Spectral Function 

The common theme in the spectral functions is that all 
dependence on the interchain coupling (and hence all im- 
portant temperature dependences in the neighborhood of 
T c ) are expressible in terms of the single coherence scale, 
A c . Moreover, it is the spectral weight of the coher- 
ent features in the spectrum, rather than their energies, 
which are strongly temperature dependent! This is very 
different from the behavior of the spectral functions near 
T c in a three dimensional BCS superconductor. 

Characteristic shapes of the single hole spectral func- 
tion above and below T c are shown in Fig. 2. Above T c , 
the single hole spectral function is a broad incoherent 
peak. Below T c , there is a coherent delta- function piece 
and a multiparticle continuum at higher energy, 

G< (fc, u) = Z{k\\ )S[w - £(k)} + G multi , (72) 

where 

£{k) = ^v 2 s k\ + Al + 2t±Z(k\i ) cos(fc_L<z) + ... . (73) 

Here kp and k± are, respectively, the components of 
the crystal momentum parallel and perpendicular to the 
chain direction. 
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FIG. 2. The temperature evolution of the spectral func- 
tion. The dashed line depicts G < (kF,ui) at temperature 
T = A s /3 > T c , and is calculated using the parameters 
7c = 0.3, K s — 1/2 and v s /v c — 0.2. The solid line rep- 
resents the spectral function at zero temperature. A coherent 
<5-function peak onsets near T c at energy Ao = A 3 + A c (0)/2. 
Here we assume A s /A c (0) = 5. The multi-particle piece 
starts at a threshold 2A C (0) away from the coherent peak. 
The exact shape of the incoherent piece at T = is not cal- 
culated in the present work and is meant to be schematic. 



The energy gap for the coherent peak is 

A (T) = A s + \a c {T) , (74) 
and its spectral weight is given by 

Z(k) ~ [A C (T)] 2 ^ . (75) 

Thus, Z(k) (and with it the transverse bandwidth) is 
the most strongly temperature dependent feature of the 
spectral function. 

The multiparticle incoherent piece G multl starts at a 
threshold energy £(k) + 2A C (T). This is the origin of 
the gap between the coherent peak and the incoherent 
shoulder in Fig. 2. Various forms of damping, including 
phase fluctuations transverse to the stripes, will broaden 
this structure, leading to a peak-dip-shoulder form of the 
spectral function. However, the distance from the coher- 
ent peak to the dip should be proportional to A C (T) and 
hence, at T = 0, to T c . 

3. The Spin Response Function 

The spin response function is entirely a multiparticle 
continuum; even below T c , we find that any spin-1 mode 
is unstable to decay into two spin 1/2 "quasiparticles" . 
However, at low temperature, we find that there is a 
spin-1 resonant state with an exponentially long lifetime 
near the threshold energy 2A S + A c = 2A , with momen- 
tum 2k p, where kp is the Fermi momentum on a stripe. 
Even here, because the barrier to decay is quantitatively 
small compared to T c , we expect that no sharp resonant 
state will appear in the spectrum in the neighborhood of 
T c . Rather, it will appear as the temperature falls below 
T « 0.2A C (0) w 0.4T C . 

B. Implications of Two Scales 

The existence of two scales in the superconducting 
state appears in different experiments in fairly obvious 
ways: 

1) Since an electron has spin and charge, the gap mea- 
sured in single particle spectroscopies, such as ARPES 
or tunnelling, is A = A s + (1/2)A C (T). (See Eq. @.) 
Manifestly, this gap scale decreases slightly with increas- 
ing temperature, but remains large, roughly A s , above 
T c . The gap scale A s is unrelated to T c , and moreover 
A S (T = 0) ^> 2T C , which physically is the statement that 
the onset of phase coherence, not pairing, is what deter- 
mines T c . Consequently, the zero-temperature superfluid 
density is a better predictorQ'Q of T c than Ao(T = 0). 
(See Eq. ([l9]) and subsequent discussion.) Similarly, 
pure spin probes, such as NMR or neutron scattering, 
see a gap which is approximately A s per spin 1/2. (See 
Eqs. ©.) 
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2) Experiments involving singlet pairs .-of electrons, 
such as Andreev tunnelling, could exhibitH an energy 
scale A c ; a scale, moreover, which vanishes at (or near) 
T c , and is related in magnitude to T c in a more or less fa- 
miliar manner, A c (0)/2 ~ T c . More complicated spectro- 
scopies, such as SIS tunnelling {e.g. tunnelling across a 
break junction) should reveal gap-like features with both 
energy scales, A s and A c . 

3) The existence of two correlation lengths implies 
that different measurements will find the order parameter 
magnitude depressed over distinct distances: If an impu- 
rity destroys the superconducting gap locally, the single- 
particle density of states, as determined, for instance, 
with a scanning tunnelling microscope, will basically re- 
cover over a length scale £ s (although, subtle effects will 
persist out to a scale £ c ). By contrast, the magnetic field 
strength near the core of a vortex, which otherwise would 
diverge logarithmically at short distances, is reduced in- 
side a "core radius" due to the fact that the superfluid 
density is depressed, (i.e. there is a lower current density 
per unit phase gradient). Since this latter effect involves 
only charge motion, the vortex core radius is of order £ c . 
This ^magnetic" core radius is measured, in principle, in 

4) The superconducting state reflects the non-Fermi 
liquid character of the normal state in many ways, but it 
has a complex scalar order parameter as in a conventional 
(BCS) superconducting state. This means that we might 
expect well-defiipd elementary excitations with the quan- 
tum numbersEjO of the electron quasiparticle, as indeed 
we have found. However, in a conventional superconduc- 
tor, the quasiparticle energy is shifted by the opening of 
the gap, and the lifetimes of all elementary excitations 
(as observed e.g. in ultrasonic attenuation) are strongly 
temperature dependent below T c . In the present case, it 
is the spectral weight associated with the elementary ex- 
citations which is strongly temperature dependent, not 
the lifetime or the energy. Moreover, even as T — > 0, the 
quasiparticle weight remains small, in proportion to a 
positive power of the distance from the quantum critical 
point. See Eq. ©. 



C. Two Scales in the High Temperature 
Superconductors 

It has been noted,@@ in the so-called "Yamada plot," 
that T c in underdoped high temperature superconductors 
is proportional to the observed incommensurabilityc3'E3 
in the low energy spin structure factor. Th e . m agnetic 
incommensurability, w/d, is inversely relatecEK3 to the 
mean separation between charge "stripes," d. Thus, the 
Yamada plot implies that T c is inversely proportional to 
the mean spacing between stripes; As the stripes become 
more separated, and the electronic structure becomes 
more one dimensional, T c -> 0. This observation strongly 
supports the idea that the anomalous electronic proper- 
ties of these materials reflect the properties of nearby 



phases of the 1DEG. Indeed, many of the spectral fea- 
tures listed above have been observed, with various levels 
of confidence, in experiments on the high temperature su- 
perconductors: 

1) The best single particle spectra (ARPES and tun- 
nelling) exist for Bi 2 Sr2CaCu208 because it cleaves eas- 
ily. For underdoped and optimally doped materials, 
the single particle gap as measured by tunnelling and 
ARPESO is found to be large: A > 35meV, in the "flat- 
band" region near the M, or (w,0) and (0, n), points of 
the Brillouin zone. (Ao/2T c lies in the range 2 to 6. 
The M region is where the maximum of a d-wave_£u- 
perconducting gap is expected.) The gap persists^ in 
some form or other to temperatures well above T c . More- 
over, Ao increases with underdoping while T c decreases. 
This gap is quite clearly a superconducting gan in that it 
has (at low T) the characteristic d-wave forroKJ expected 
of the superconducting gap, and it evolvesE-3 smoothly 
with overdoping into a gap of only slightly smaller mag- 
nitude which opens, in a more conventional manner, in 
the neighborhood of T c . (We focus on the gap near the. 



M point, especia 
and experimental 



because there are both theoreticaEI 
reasons to think that the "flat-band" 
region is associated with states in "stripes" or "fluctuat- 
ing stripes.") Moreover, as discovered first by Uemura 
and coworkers,!! T c is roughly proportional to the zero 
temperature superfluid density for underdoped materi- 
als, consistent with the notiontj that it is phase ordering, 
not pairing, whkh determines T c . 

2) DeutschertHl has argued that the gap scale de- 
termined by low temperature Andreev tunnelling spec- 
troscopy is considerably smaller than that determined 
from single particle tunnelling measurements in under- 
doped materials, while the two gap scales approach each 
other in overdoped materials. This issue is well worth 
revisiting in more detail. The single particle gap scale 
is strongly apparentED in SIS tunnelling spectra - we do 
not know of any convincing analysis which reveals the 
smaller charge gap scale in such experiments. 

3) The vortex core radius has hepu measured with both 
scanning tunnelling microscopyMO (STM) and /xSR.EJ 
The /iSR study measures the magnetic field distribution 
in the material, and infers the core radius from the high- 
field cutoff of the distribution. For large applied fields 
(B ~ 6T), both methods are in rough agreement that 
the core radius is about 15A. However, the core radius 
deduced from the /xSR measurements is strongly field de- 
pendent, so that at low fields (B ~ 0.5T) it yields a cc*a 
radius around 120A. By contrast, preliminary evidenceMI 
from STM experiments suggests that the core radius 
measured by that method is not strongly field dependent, 
so that, in low fields, the results of the two methods differ 
by almost an order of magnitude. However, there appeax. 
to be differences in the STM results of different groups.. l£ 
Certainly, the core radius inferred from STM studieala-EH 
of the gap suppression in the vicinity of an impurity at 
zero magnetic field are suggestive of a rather short co- 
herence length. While the experimental results are, by 
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no means, definitive, we would tentatively like to explain 
the discrepancy between the STM and fiSR results at low 
field as evidence of the existence of two coherence lengths 
in the superconducting state. 

4) It has been realized for a long time that there are no 
sharp quasiparticle features in the ARPES spectrum near 
the superconducting gap maximum (near the M point of 
the BrUJLouin zone) in the normal state, and it has been 
arguedEJ that they disappear due to a lifetime catastro- 
phe which occurs as the temperature is raised above T c . 
RecentcrEl high-resolution ARPES measurements in opti- 
mally doped P^S^CaC^Os have revealed a new picture 
of the emergence of these peaks. Within experimental 
resolution, neither the energy nor the width of the peak 
changes as the temperature is raised from well below T c 
to slightly above T c ; rather, it is the intensity of the peak 
that is strongly temperature dependent in the neighbor- 
hood of T c . The intensity vanishes slightly above T c , 
without any apparent change in the shape of the peak 
itself. Indeed, the sharp temperature dependence of this 
intensity in the neighborhood of T c is consistent with its 
being proportional to a fractional power of the (local) 
superfiuid density. (See Eq. (|75|).) Additional evidence 
for this comes from an old observation of Harris et alSB 
that, as a function of undcrdoping, the weight in the 
peak at low temperatures decreases with decreasing su- 
perfiuid density. Moreover, Shen and BalatskyEj have 
argued that a small dispersion of the ARPES peak in the 
direction perpendicular to the putative stripe direction 
scales more or less with T c , consistent with our Eq. (61). 
The distance between the coherent peak_and the dip fea- 
ture in ARPES curves near the M pointEJ decreases with 
underdoping, consistent with the zero temperature dis- 
tance being proportional to T c . 

A similar temperature evolution has been observed-ifor 
the so called "resonant peak" in neutron scatterings in 
YBa 2 Cu307_«5 and Bi^Si^CaC^Og (although no such 
feature has been seenE3 in L^-zSr^CuC^). We would 
like to identify this phenomenon, as well, with a dimen- 
sional crossover of the sort discussed here. However, the 
spin resonance we have found in the present model is 
clearly not directly related to the observed resonant peak. 
In particular, all features we have found are peaked at a 
momentum 2kp, while the resonant peak is centered on 
the antiferromagnetic wave- vector (tt,tt). Moreover, the 
peak we have found disappears through a lifetime catas- 
trophe well below T c , while the resonant peak is sharp 
immediately below T c . Clearly, at the very least, to have 
a theory with anything more than a very rough caricature 
of this observed magnetic behavior, we need to expand 
the considered modeEl to include the effects of the anti- 
ferromagnetic "strips" between the "stripes." 



D. Further Implications for High Temperature 
Superconductors 

Finally, we end with a few additional observations con- 
cerning insights into the behavior of the high temperature 
superconductors that can be obtained from the analysis 
of this paper: 

In the superconducting state of a structurally quasi- 
1D superconductor, we have found there are two emer- 
gent length scales, £ s and { c . If the quasi- ID electronic 
structure is self organized, as it is in the high temper- 
ature superconductors, there p*e potentially two addi- 
tional emergent length scales:E-3 the mean spacing be- 
tween stripes, d, and the persistence length of the stripes, 
£,stripe- d can be determined directly ..from the charge 
incommensurabilitjO'Ej (or indirectlyoCj from the spin) 
structure factor. Stripe is much harder to determine ex- 
perimentally, although it is bounded below by the cor- 
relation length of the magnetic order So long as 
{stripe is the longest length scale in the problem, i.e. 
so long as {stripe £o it is possible to assume, as we 
have here, that the superconducting properties of the 
system are quasi one-dimensional. Where this inequal- 
ity is violated, the correct theory of the superconduct- 
ing state needs to be .significantly modified. As was 
pointed out previously,^ so long as the weaker condi- 
tion, {stripe 3> { s , is satisfied, it is possible to have a one 
dimensional theory of spin gap formation. At present 
experiments are unclear about the range of doping and 
materials for which either of these inequalities is satis- 
fied, which is the most important source of uncertainty 
in the application of these ideas to the high temperature 
superconductors. Certainly, with sufficient overdoping, 
the stripes loose their integrity and the application of 
these ideas becomes suspect. 

To get a feeling for magnitudes, we can make rough 
quantitative estimates of the remaining length scales 
from well-established experimental data in the high tem- 
perature superconductors, although numbers vary from 
material to material, and as a function of doping con- 
centration, x. The spin velocity in the undoped anti- 
ferromagnet is around v s rs 0.8eV-A, and the supercon- 
ducting gap is A ~ 35meV, so £ s - 20A. The charge 
coherence length { c = { s (v c /v s )(A s / A c ), so if we esti- 
mate v c /v s ~ t/J w 2 — 3 and A c w 2T C « 16meV, 
we find that characteristically { c ~ 100 — 150A._(|This is 
in good agreement with the fiSR measurements of the 
vortex core radius cited above. ) The spacing between 
stripes is in the range of four or more lattice constants, 
d ~ 16A. 

A crossover magnetic field, which can be identified as 
a mean fieldEa B C 2 , can be estimated as the field at which 
there is one vortex per coherence length { c between each 
pair of neighboring stripes; this leads to an estimate 

B c2 ~ <f>o/{ c d , (76) 

where (j)Q — hc/2e is the superconducting flux quan- 
tum. While B C 2 estimated in this fashion is quite large 
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(cj) /£ c d = 80T for d = 16Aand £ c = 100A) it is small 
compared to the characteristic magnetic field, 

B s ~ 2fo/t.w , (77) 

at which orbital effects lead to the destruction of the spin 
gap. Here, w is the "width" of a stripe - i.e. the width 
of the Id region involved in spin gap fownation. In the 
"spin gap proximity effect" mechanismlia proposed pre- 
viously this would imply that w is one to two times the 
crystalline lattice constant. The extremely large value 
of B s rationalizes the lack of any observable reduction 
of the spin gap temperature in the recent NMR exper- 
iments of Gorny et alSB up to fields as high as 12T in 
YBa2Cu307_5. (See, also, the discussion in Ref. |82] .) 

One important difference between a stripe phase and 
the array of chains studied here, is that in the stripe phase 
there are additional electronic degrees of freedom which 
live in the antiferromagnetic strips between the stripes. 
The two-component nature of the electronic structure of 
doped antiferromagnets,E3 is characteristic of the micro- 
phase separation physics that gives rise to this state. Of 
course, the antiferromagnetic strips are themselves quasi 
one-dimensional magnets,0 so that any magnetic order- 
ing must be viewed, in similar spirit to that considered 
here, as resulting from a dimensional crossover. Indeed, 
it is certainly the spins in tie insulating strips that make 
the dominant contributionEJ to the "resonant peak" ob- 
served in neutron scattering. A detailed theory of this 
peak is beyond the scope of the present model, but is 
embodied in the spin gap proximity effect .113 

However, we have found a new neutron scattering reso- 
nance for a quasi one-dimensional superconductor. While 
the dimensional crossover causes no bound state in the 
spin-1 excitations, we find a resonant state of two spin-^ 
quasiparticles appearing below T m 0AT c . The mode ap- 
pears at an energy 2 A,, + A c = 2Ao, or twice the single 
particle gap as measured by ARPES or tunneling, and 
at momentum 2fcp, where fcp defines the Fermi surface 
associated with a stripe. Since this is a four soliton res- 
onance, it may be qualitatively sensitive to deviations 
from the limit A c < < A s , so that the resonance is likely 
to be most well defined in the underdoped region where 
T C «A . 

Finally, we remark that the ARPES spectrum along 
the symmetry direction from (0,0) to i.e. along 

the ray which is expected to pass through the node of a 
d-wave gap function, is very different in character from 
that in the M that we have discussed. In clean saniales of 
optimally doped Bi 2 Sr 2 CaCu208, there is a peakOO in 
the spectral function both above and below T c , and the 
peak reaches the Fermi surface at a well defined "nodal 
point", k n = (0.447T, 0.447r). This peak does,not exhibit 
the characteristics of a quasiparticle peakjla in that its 
width is always larger than its energy; indeed, it seems to 
exhibit quantum critical behavior reminiscent of a LuU 
tinger Liquid. Moreover, there is no qualitative changeO 
in the temperature evolution of this peak as the tem- 
perature is lowered from two or three times T c down to 



temperatures as low as at least 1/2 T c ; the character of 
the nodal excitation seems to be remarkably insensitive 
to the onset of superconductivity. By contrast—in opti- 
mally doped La2_ a; Sr 2 ;Cu04there is apparentlyO no ob- 
servable peak along the nodal direction, and indeed little 
or no spectral weight within about 0.5 eV of the Fermi 
energy. Indeed, recent neutron scattering studiesO of 
the low energy magnetic scattering in the neighborhood 
of 2k n have revealed the existence of a clean spin gap at 
low temperatures, which is apparently inconsistent with 
the existence of any gapless nodal quasiparticle excita- 
tions. 

It is clear that whatever spectral response is observed 
near k n is not associated with the vertical and hori- 
zontal stripes studied here, because a stripe wave vec- 
tor does not span the "Fermi surface" along this di- 
rection. It could be associated with diagonal stripes, 
which have-,been observed recently in various insulating 
materials,! 76 ^ 7 ! in which case the observed quantum criti- 
cal behavior might truly be that of a Luttinger liquid. An 
alternative picture is backflow associated with holes that 
have not condensed into vertical or horizontal stripes. 
Both explanations are conceivable as there are strong rea- 
sons to expect that the orienting potential, which locks 
the stripes along a particular (vertical or horizontal) crys- 
tallographic direction to be stronger in La2-^Sr x Cu04 
than in Bi2Sr2CaCu20s. However, other sources of quan- 
tum critical behavior are certainly possibleH We will de- 
fer further discussion of these classes of excitations to a 
future study. 
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APPENDIX A: THE EFFECT OF INTERCHAIN 
SINGLE PARTICLE HOPPING 

Until now, we have ignored single particle hopping be- 
tween chains. This is because, especially in the presence 
of a spin gap, it is irrelevant in the renormalization group 
sense. However, in the superconducting state, we expect 
the quasiparticles to be able to propagate coherently be- 
tween chains. Because these terms are irrelevant, their 
effects on the spectrum can be computed in ordinary de- 
generate perturbation theory. It is easy to see that to 
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first order in the interchain hopping, the quasiparticle 
energy is 

£{k) = y^fAf + A§ + Z(Jfe||)eW(*±) + 0(e { ^f , (Al) 

where e^- ± \k±) — 2tj_ cos(fcj_a) is the interchain contri- 
bution to the quasiparticle dispersion, and k\\ +kp and fcj_ 
are, respectively, the components of the crystal momen- 
tum parallel and perpendicular to the chain direction. 

This is highly reminiscent of the spectrum we would 
have obtained were we to compute the spectrum of a 
quasi one-dimensional superconductor using BCS mean 
field theory 



g(BCS) 



[v F k + e^)(k ± )} 2 + A 2 (A2) 
= £[ BCS) + [v s k/£[ BCs y^(k^) + 0( e W) 2 , 



with the differences that the Fermi velocity is replaced by 
the (slower) spin velocity, the superconducting gap is the 
sum of the (single-chain) spin gap and the (interchain) 
charge gap, and the interchain bandwidth is reduced by 
the quasiparticle weight factor Z. 



APPENDIX B: MACROSCOPIC SUPERFLUID 
DENSITY 

In this appendix, we compute the macroscopic phase 
stiffness (superfluid density) tensor i^ a b[«;] in two dimen- 
sions (a = x, y) given a microscopic distribution of the (in 
general anisotropic) local phase stiffness tensor, K ab {r)- 
We include the derivation here for pedagogical purposes, 
although the results exist elsewhere in the literatureE-i 

k determines the relation between the local current 
density, j(f) and the gradient of the phase according to 



Ja{r) 



n ab d b 8(f) 



(Bl) 



From the equation of continuity, it follows that V • j = 0, 
so we can express j in terms of a potential, j a (r) = 
£ a bdb4>{r), so that 



£abdb(j>{r) = n ab d b 8(r) . 



(B2) 



To compute K xx in a rectangular geometry, this equation 
is to be solved subject to the boundary conditions that 
8 = for x = and 8 = A8 for x = L x (independent of 
y) and (from the condition that no current can flow out 
of the sample in the y direction) = for y = and 
<f> = A<1> for y = L y . For a given distribution of k, we 
solve this equation for given A8 to determine A<p, from 
which we determine K according to 



K, 



Acj>/A8 



(B3) 



The key obser vati on is the same potential and phase 
that satisfy Eq. (B3), also satisfy the dual equation 



where 



Therefore 



£acKj(r)e db 



K xx [n]K y 



= 1 



(B4) 



(B5) 



(B6) 



We can apply this general result to the problem of in- 
terest here. Consider the case of a square geometry in 
which, because of some assumed domain structure, the 
system is macroscopically isotropic (R = K xx = K yy ) 
despite the existence of microscopic anisotropy in each 
"stripe" domain. It follows that 

k(kx-> /cm ) k(1/k», = 1 • (B7) 

It follows that k(kj^,K||) = ^k±k\\ . Other solutions to 



Eq. (B7) exist, but are not homogeneous functions. 



APPENDIX C: THE EFFECTIVE POTENTIAL, 

V(X) 

To compute the effective potential which appears 
in Eq. (|62|), we considec_.the discrete version of the 
refermionized HamiltonianEH, 

H = -]>> + (-!)" A(n)/2Mc„ +1 +H.C.] , (CI) 



where v r = 2tn and 



A(n) = A c sign(i? 2 - An 2 ) , 



(C2) 



corresponding to a pair of solitons separated by a dis- 
tance R. (We have set the lattice constant equal to 1.) 

We compute the ground state energy on a system of 2N 
sites by computing the single particle eigenvalues, and 
then summing over the lowest lying N — 1 of them to get 
the total energy as a function of R. This is precisely the 
program carried out previously to study iiariaus solitonic 
states in the SSH model of polyacetylene.EJ~E£l With open 
boundary conditions, the Hamiltonian matrix is tridiag- 
onal, and so particularly simple to study numerically for 
large system sizes. We have carried out this program nu- 
merically for system sizes up to 2N — 3000, and for A c = 
0.2, 0.1, 0.05, and 0.02; the continuum limit is obtained 
when A c — > 0. Even for the smallest values of A c , we 
find no significant finite size effects at these large system 
sizes. The results for v(x) computed in this way are sum- 
marized in Fig. 1. The fact that the asymptotic value of 
v is always slightly negative is a reflection of the fact that 
in the limit of small soliton width (which is equal to 1, the 
lattice constant, in the present calculation), the soliton 
creation energy is a very strongly varying function of the 
width— as found previously by Takayama, Lin-Liu, and 
MakiL3, and only approaches its true asymptotic limit 
A c /2, when l/£ c is extremely small. 
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